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We present new experimental data on directed flow in colli- 
sions of Au+Au, Xe+Csl and Ni+Ni at incident energies from 
90 to 400A MeV. We study the centrality and system depen- 
dence of integral and differential directed flow for particles 
selected according to charge. All the features of the experi- 
mental data are compared with Isospin Quantum Molecular 
Dynamics (IQMD) model calculations in an attempt to ex- 
tract information about the nuclear matter equation of state 
(EoS). We show that the combination of rapidity and trans- 
verse momentum analysis of directed flow allow to disentangle 
various parametrizations in the model. At 400^4 MeV, a soft 
EoS with momentum dependent interactions is best suited 
to explain the experimental data in Au+Au and Xe+Csl, 
but in case of Ni+Ni the model underpredicts flow for any 
EoS. At 90 A MeV incident beam energy, none of the IQMD 
parametrizations studied here is able to consistently explain 
the experimental data. 

PACS: 25.70.Lm, 21.65.+f, 25.75.Ld 



I. INTRODUCTION 

The study of collective flow in relativistic heavy-ion 
collisions has been an intense field of research for the 
past twenty years (see Refs. [1,2] for recent reviews). The 
ultimate motivation for the whole endeavour has been 
the extraction of the equation of state (EoS) of nuclear 
matter (see Ref. [3] for an early account and Ref. [4] for 
more recent ones) . Moreover, the study of highly complex 
(quantum) many-body dynamics of heavy-ion collisions 
is in itself a challenging task. 

The (in-plane) directed (or sideward) flow was pre- 



dicted for semi-central heavy-ion collisions on the ba- 
sis of fluid dynamical calculations [5] and observed in 
experiments soon after [6,7]. The study of the average 
in-plane transverse momentum, (p x ), as a function of ra- 
pidity, y, provides an easy and intuitive way of quan- 
tizing the directed flow [8]. For beam energy range up 
to a few GeV per nucleon, the experimental [8-24] and 
theoretical [25-43] studies of this phenomenon have pro- 
vided important understanding of its features. Exper- 
imentally, the dependences of directed flow on incident 
energy [9,16,19,18,22], centrality [9,15,17,19-22], particle 
type [10,16,17,19-22], system size [9,18,22] and isospin 
[21] have been determined. These main features have 
been quantitatively reproduced by microscopic transport 
models of either Quantum Molecular Dynamics (QMD) 
[29,33,15,16,20,41,24] or Boltzmann-Uehling-Uhlenbeck 
(BUU) [11,30,31,36,38,17,40,42,43,24] type. Despite all 
these efforts, a definite conclusion on the EoS has not 
yet been achieved. As pointed out early on [27,31,32], 
the momentum dependent interactions (MDI) play a cru- 
cial role in the determination of the EoS. The directed 
flow is influenced in addition by the (in-medium) nucleon- 
nucleon cross section (a nn ) [36,37,44]. Moreover, consis- 
tency is needed in deriving EoS together with both MDI 
and a nn [35,43]. The importance to combine various ob- 
servables, each sensitive (ideally) to one particular aspect 
of the parameters that influence EoS was emphasized re- 
cently [4]. 

The results quoted above have been obtained from the 
analysis of integrated directed flow. The first analysis of 
differential directed flow (DDF) was introduced by Pan 
and Daniclewicz [36] , who studied the transverse momen- 
tum (p t ) dependence of the first order Fourier coefficient, 
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Vi = (cos((j))), where <f> is the angle with respect to the re- 
action plane. The DDF was studied around the balance 
energy (E{,ai , which is the energy of disappearance of flow 
[14,16,23]) by Li and Sustich [44], who unraveled its inter- 
esting patterns. They also pointed out the marked sen- 
sitivity of DDF to both EoS and nucleon-nucleon cross 
section (cr nn ). At AGS energies the DDF was studied 
both experimentally [45] and theoretically [46,47]. Re- 
cently, we have completed the first experimental analysis 
of DDF for Au+Au collisions at incident energies from 90 
to 400 A MeV [48] . We have found interesting patterns of 
the differential flow, evolving as a function of incident en- 
ergy, particle type and rapidity. In particular, the study 
of high-p t particles is important because, as proposed in 
Ref . [49] , they are good messengers from the high density 
state of the collision. The DDF could additionally pro- 
vide snapshots of the flow development during the time 
of the collision. 

In this paper we present new experimental data on di- 
rected flow in collisions of Au+Au, Xe+Csl and Ni+Ni. 
The complete coverage of the FOPI detector makes pos- 
sible precision studies of flow, refining our earlier studies 
done with Phase I data [15,19,20,22]. Following our re- 
cent exploration of the energy range of 90 to 400A MeV 
for Au+Au [48], we focus here on the centrality (for 
Au+Au) and on system size dependence of directed flow. 
This analysis has been performed for the incident ener- 
gies of 250 and 400 A MeV, for particles with Z=l and 
Z=2. To avoid overloading, only a selection of results is 
included in the body of the paper. In the Appendix we 
provide additional figures to complete the data set. After 
describing the detector, the method of analysis, and the 
corrections applied to data, we study the centrality and 
system dependence of directed flow over the complete for- 
ward rapidity range, both in terms of (pt) integrated and 
in a differential way. All the features of the experimental 
data are then compared with IQMD model calculations 
for the incident energies of 90 and 4QQA MeV. 

II. SET-UP AND DATA ANALYSIS 

The data have been measured with a wide phase-space 
coverage using the FOPI detector [50] at GSI Darmstadt. 
The reaction products are identified by charge (Z) in the 
forward Plastic Wall (PW) at 1.2° < d iah < 30° using 
time-of-fiight (ToF) and specific energy loss. In the Cen- 
tral Drift Chamber (CDC), covering 34° < 0i ab < 145°, 
the particle identification is on mass (mass number A), 
obtained using magnetic rigidity and the energy loss. For 
PW the Z resolution is 0.13 charge units for Z~l and 
0.14 for Z=2, while for CDC the mass resolution varies 
from 0.20 to 0.53 mass units for A~l to A=4. The con- 
tamination of Z=l in the Z=2 sample varies from 6% to 
10% (from Ni+Ni at 250A MeV to Au+Au at 400,4 MeV) 
for the PW and is up to 20% for the CDC (where it 
is the contamination of A— 1,2 and 3 in the A=4 sam- 



ple) . The PW measures the velocity of particles via ToF 
with an average resolution of 150 ps. For the CDC, the 
relative momentum resolution cr Pt /p t varies from 4% for 
p t < 0.5 GeV/c to about 12% for p t =2 GcV/c. For more 
details on the detector configuration for this experiment 
see Ref. [51]. 




FIG. 1. FOPI detector acceptance: phase-space distribu- 
tion for Z=l particles measured in the centrality bin M4 of 
the reaction Au+Au at 250A MeV. The intensity contours are 
spaced logarithmically. The thicker lines mark the geometri- 
cal acceptance of different subdetectors. 

The phase-space coverage of the FOPI detector is pre- 
sented in Fig. 1 for particles with Z=l (for PW) and 
A=l,2,3 (for CDC), measured in semi-central collisions 
Au+Au at incident energy of 250A MeV. To compare 
different incident energies and particle species, we use 
normalized center-of-mass (cm.) transverse momentum 
(per nucleon) and rapidity, defined as 

p{ 0) = ( Pt /A)/(p^-/A P ), y(°) = (y/y P y m -, 

where the subscript P denotes the projectile. For the 
PW coverage, shadows around 6i a b = 7° and 19° are 
visible, arising from subdetector borders and frames. 

For the centrality selection we used the charged parti- 
cle multiplicities, classified into five bins, Ml to M5. The 
variable Erat = £a,i/ JV E« , (the sums run over the 
transverse and longitudinal cm. kinetic energy compo- 
nents of all the products detected in an event) has been 
additionally used for a better selection of the most cen- 
tral collisions (M5 centrality bin) . The geometric impact 
parameters interval for the centrality bins M3, M4, and 
M5 for Au+Au system at 400,4 MeV studied here are 
presented in Table I. 

The impact parameter intervals corresponding to the 
three investigated systems at the incident energy of 
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TABLE I. The geometric impact parameters intervals 
Ab geo and the correction factors for the reaction plane resolu- 
tion, l/(cos Ac/)}, for three centrality bins of Au+Au collisions 
at the incident energy of 400A MeV. 



Centrality bin 


M3 


M4 


M5 


Ab geo (fm) 


6.1-7.6 


1.9-6.1 


0-1.9 


l/{cos A4>) 


1.05 


1.04 


1.17 



TABLE II. The geometric impact parameters intervals 
Abg eo , the reduced impact parameters (b ge o) /f)^" 1 and 
the correction factors for the reaction plane resolution, 
l/(cosA0), for the three systems at the incident energy of 
250A MeV, M4 centrality bin. 



System 


Au+Au 


Xc+Csl 


Ni+Ni 


Abg e(> (fm) 


1.9-6.1 


1.7-4.8 


1.5-3.4 


lb \lh max 


0.31 


0.29 


0.27 


l/(cos A(j>) 


1.05 


1.09 


1.27 



250A MeV, M4 centrality bin, are presented in Table II 
along with the reduced impact parameters (b geo ) /b™ a Q x . 
is the maximum geometrical impact parameter, cal- 



hrnax 
geo 

culated as: b™ x 



1.2(Ap /3 + A]/ 3 ) (in fm). Applying 
the same recipe for the centrality selection for different 
systems the reduced impact parameter is similar, as seen 
in Table II. 

The largest data sample was acquired with what we 
call "Medium bias" trigger, which accepts events roughly 
corresponding to centrality bins M3, M4 and M5 for 
Au+Au collisions. For Xe and Ni systems this trigger 
selection amounts to a bias for the M3 centrality. This is 
the reason why M3 is not included in the present paper 
for these systems. We have collected "Minimum bias" 
data for all systems, but the statistics is far smaller, thus 
not allowing the type of analysis done in this paper. 



A. The reaction plane determination and the 
correction for its resolution 

The reaction plane has been reconstructed event-by- 
event using the transverse momentum method [8]. All 
charged particles detected in an event have been used, 
except a window around midrapidity (\y^\ < 0.3) to 
improve the resolution. The particle-of-interest has been 
excluded to prevent autocorrelations. 

The correction of the extracted values due to the recon- 
structed reaction plane fluctuations has been done using 
the recipe of Ollitrault [52]. The resolution of the reac- 
tion plane azimuth, A<j>, can be extracted by randomly 
dividing each event in two subevents and calculating for 
each one the reaction plane orientation, <!>i and $2 [8,52]. 
From the resolution, quantified as (cos($i — $2)), the 
correction factors, 1 / (cos A^) can be calculated [52] (see 
Rcf. [51] for more technical details). For the experimen- 
tal data the correction factors for the centrality bins M3, 



M4, and M5 for Au+Au system at 400A MeV are pre- 
sented in Table I. The values for the three investigated 
systems at the incident energy of 25QA MeV, M4 central- 
ity bin, are presented in Table II. 



IQMD HM Au+Au E=400 AMeV 
• 1/<COS(A(]))> 
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FIG. 2. Upper panel: the resolution of the reconstructed 
reaction plane (squares) and the corresponding correction fac- 
tors (dots). Lower panel: v\ values for the true (continuous 
line) and reconstructed and corrected (dashed line) reaction 
plane. IQMD HM events were used for these studies. 

The accuracy of the reaction plane correction proce- 
dure was checked using Isospin Quantum Molecular Dy- 
namics (IQMD) [33] events analyzed in the same way 
as the data. The results are shown in Fig. 2 for IQMD 
events at 400A MeV. The upper panel presents the res- 
olution and the correction factor as a function of the 
impact parameter, b. Their dependence on centrality 
reflects mainly the dependence of the strength of the 
directed flow (see lower panel), but also a finite num- 
ber effect [8], evident towards peripheral collisions. The 
lower panel of Fig. 2 presents the centrality dependence 
of integrated v\ values for Z=l particles (for the forward 
hemisphere), derived from IQMD events in two cases: i) 
with respect to the true reaction plane (known in the 
model) and ii) with respect to the reconstructed reaction 
plane and using the correction according to [52] (correc- 
tion factors of the upper panel of Fig. 2). The agreement 
between the two cases is perfect, down to most peripheral 
collisions, where correction factors up to 2 are necessary. 

Alternative methods of flow analysis have been pro- 
posed recently [53]. However, because, for our energy 
domain, the flow of nucleons is at its maximum and pro- 
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duced particles are very rare, the impact of these refined 
methods is expected to be minor. 



B. The influence of FOPI detector on the flow 
measurements 

As seen in Fig. 1, the complete phase space coverage of 
the FOPI detector (in its Phase II) is hampered by one 
empty region, corresponding to polar angles 9i a b from 
30° to 34° in the laboratory frame. Additional detector 
shadows around 7° and 19° are present too. We have 
studied the effect of the FOPI acceptance using IQMD 
transport model [33]. The IQMD events were analyzed 
in the same way as the experimental data. 
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FIG. 3. The effect of the FOPI detector filter on vi val- 
ues for Z=l and Z=2 particles for the incident energy of 
400.4 MeV, M4 centrality bin. IQMD SM events were used 
for this comparison. 

The results are presented in Fig. 3, where the inte- 
grated V\ values as a function of rapidity for Z=l and 
Z=2 particles are shown for the ideal case of total cov- 
erage (full lines) and when the FOPI filter is employed 
(dashed lines). IQMD SM events for the incident en- 
ergy of 400A MeV, M4 centrality bin were used for this 
comparison. With the present configuration of the de- 
tector, the measured V\ values are very close to the ideal 
case. We note that, for our published directed flow data 
[15,19,20,22], although the effect of the FOPI acceptance 
on the directed flow results was quite small [15,54], its 
magnitude was comparable to the difference between soft 
and hard EoS. 

An important ingredient in the present analysis is the 
correction for distortions due to multiple hit losses. As an 
example, in case of PW, despite its good granularity (512 
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FIG. 4. The effect of the FOPI detector on vi values for 
Z=l particles, for the incident energy of 400^1 MeV, M4 cen- 
trality bin. IQMD events were used for this comparison. 
Three cases are compared: the simple geometric FOPI fil- 
ter (full line) and GEANT simulation without (dotted line) 
and with (dashed line) the multiple hit correction. 



independent modules [50]), average multiple hit proba- 
bilities of up to about 9% at 4QQA MeV are registered 
for the multiplicity bin M4. Because of the directed flow, 
the average number of particles detected over the full PW 
subdetector is up to 2 times higher in-plane than out of 
the reaction plane. As a consequence, the losses due to 
multiple hit are strongly correlated with the directed flow 
and follow its dependences on incident energy, centrality 
and system size. These losses lead to an underestimation 
of the measured directed flow and need to be taken into 
account. We developed a correction procedure based on 
the experimental data, by exploiting the DDF left-right 
symmetry with respect to midrapidity. The correction 
acts upon v\ values (namely only on average values) and 
is deduced for each system, energy and centrality sepa- 
rately. Due to the flow profile in the polar angle, it de- 
pends also on transverse momentum. It is larger for Z=2 
compared to Z—\ particles. The correction was derived 
in a window around midrapidity (I?/ - 1 < 0.1) and prop- 
agated for other rapidity windows for each pf bin along 
the lines of constant 9i a b to follow the detector segmen- 
tation. It reaches up to 12% for Au+Au at AQQA MeV 
and is almost negligible at 90A MeV. At 400A MeV it is 
up to 5% for Xe+Csl and up to 2% for Ni+Ni. 

The procedure was checked and validated using IQMD 
events passed through a complete GEANT [55] simula- 
tion of the detector. IQMD events were used for this 
study, at the incident energy of 400A MeV. The results 
are presented in Fig. 4 for Z—\ particles in M4 central- 
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ity bin. The Vi values extracted from IQMD events as 
inputs into a complete GEANT simulation of the detec- 
tor (and analyzed in exactly the same way as the data), 
without (dotted line) and with (dashed line) the multi- 
ple hit correction, are compared with the true v\ values, 
obtained from standard IQMD events (actually the same 
events used for the GEANT simulation) with only a sim- 
ple geometric FOPI filter (full line). It is obvious that 
the correction is restoring the "true" v\ values (IQMD) 
from the "measured" GEANT values. Also, the correc- 
tion used for the data is quantitatively reproduced by 
these simulations. Note that, because of larger multi- 
plicities for Z—l particles from IQMD (see Section IV), 
the correction is larger in the simulations than in the 
data. As a result of these studies, all the experimental 
data (both differential and integrated) have been cor- 
rected according to the procedure described above. 

The only source of systematic error on our measured 
v\ values could be the correction for multiple hit losses 
outlined above. However, as we have demonstrated based 
on complete GEANT simulations, this correction is well 
understood. As a result, the systematic error depends on 
incident energy, centrality, particle type and p[°\ It is 
below 5% on the differential v\. There are exceptions 
for some points, for which the systematic error arises 
from (rapidity-dependent) regions in pf^'m which detec- 
tor shadows are influencing the data. For those particu- 
lar points the systematic error is already included in the 
plots. For the integrated v\ values the error is slightly 
smaller, up to 4%, including the influence of the uncov- 
ered region of 6i a b = 30° — 34°. These values do not 
include the effect of particle misidentification. 




FIG. 5. Centrality dependence of the integrated directed 
flow as a function of rapidity for Au+Au at 400A MeV for 
Z=2 particles. The lines are joining the symbols to guide the 
eye. 



III. GENERAL FEATURES OF THE DATA 

A. Centrality dependence 

By varying the centrality of the collision one aims at 
controlling both the size of the participant fireball (and 
consequently the magnitude of the achieved compression 
and subsequent expansion) and the size of the spectator 
fragment region. While semi-central collisions could pro- 
vide information preferentially on (density dependent) 
EoS, more peripheral reactions can help in pinning down 
theMDI[4]. 

Figure 5 shows the centrality dependence of the di- 
rected flow for Au+Au at A00A MeV for Z=2 parti- 
cles. Plotted as a function of rapidity are the (pt) in- 
tegrated v\ values (upper panel) and those integrated 
values weighted by the average transverse momentum, 
ip^), for the respective rapidity bin. This weighted v\ 
quantity is proportional to the average in-plane trans- 
verse momentum (p x ) (v\ = (p x /pt)) and was chosen due 
to its convenience in applying the corrections discussed 
in Section II B. First, one can notice that the known 



behaviour of the slope at midrapidity, namely the max- 
imum for intermediate impact parameters (M4 bin), is 
evident only for the weighted v\ values (lower panel of 
Fig. 5). The asymmetries {v\ values, upper panel) are 
the same around midrapidity for M3 and M4 centrality 
bins. Second, in both observables, the most significant 
dependence on centrality is taking place in the spectator 
region (roughly > 0.5) and it is more pronounced 
for the weighted v\ values. Both the asymmetries and 
the in-plane transverse momentum reflect the influence 
of the participant and of the spectator size controlled by 
the variation of the centrality. These distributions are 
inherently a result of the superposition of collective and 
thermal contributions [22]. For a given flow magnitude, 
higher temperatures (presumably achieved for more cen- 
tral collision) would translate into a smaller effective flow. 

In Fig. 6 we show the centrality dependence of the 
differential flow for Z=2 particles for Au+Au collisions 
at 400A MeV. Three centrality bins (different panels) are 
compared for three rapidity windows (different symbols). 
The lines are polynomial fits to guide the eye. The ra- 
pidity dependence of the DDF in different centrality bins 
follows the rapidity dependence of the integrated flow 
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FIG. 6. Differential flow for three centrality bins, in three 
rapidity windows, for Z=2 particles for collisions Au+Au at 
400A MeV. The lines are polynomial fits to guide the eye. 
The arrows mark the values of the average Pt for the corre- 
sponding centrality bin. 

seen in Fig. 5: the most pronounced dependence is reg- 
istered for the most central collisions, M5. As we have 
already discussed in [48] , the shape of the DDF (a gradual 
development of a limiting value, followed by a decrease 
at high pf" 1 ) could be a result of the collision dynamics. 
Part of the high-pt particles could have been emitted at 
a pre-cquilibrium stage, therefore not reaching the max- 
imum compression stage of the reaction. However, this 
possibility seems to be ruled out by the observation that 
the high-p t particles originate preferentially from high- 
density regions of the collision [49] . 

The arrows in Fig. 6 mark the values of the aver- 
age p^for the corresponding centrality bin, according to 
the symbols. For this incident energy of 400A MeV the 
value of the projectile momentum in the cm. system is 
433 MeV/c per nucleon. Higher values of average trans- 
verse momenta are seen for more central collisions as a 
result of a stronger expansion from a bigger and more 
compressed source. The dependence of (pj ^ } on rapidity 
is different for the most central collisions (M5 bin) com- 
pared to semi-central ones, for which smaller transverse 



momenta are seen towards the projectile rapidity as a 
result of the influence of the spectator matter. 

Similar data as those presented in Fig. 5 and 6 are pre- 
sented in the Appendix for Z=l particles in Au+Au at 
400A MeV and for Z=l and Z=2 particles at 250.4 MeV 
(Fig. 19 to Fig. 24). 

B. System size dependence 

As for the centrality variation, by varying the system 
size one aims to control the size of both the participant 
and the spectator. However, the question whether trans- 
parency plays a role in case of lighter systems, needs to 
be addressed simultaneously, as it results in a decrease 
of the achieved compression. In addition, the surface (or 
surface-to-volume ratio) can play an important role. It 
was suggested that the system size dependence of the 
directed flow could give insights about a nn [32]. 




FIG. 7. Phase space distributions d 2 N/dpf ) dy (0) of 
Z=2 particles for three systems at the incident energy of 
25CL4 MeV, M4 centrality. 

Figure 7 presents the phase space distribution 
d 2 N/dp^dy(°) of Z=2 particles for the three systems 
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at the incident energy of 250A MeV, centrality bin M4. 
From Au+Au to Ni+Ni system, the phase space popula- 
tion becomes more and more focused, both in transverse 
momentum and rapidity. This is an indication of the de- 
crease of stopping for lighter systems. Maximum density 
reached in the fireball depends on the system size [34], 
presumably as an effect of different stopping. On the 
other hand, due to the sizes of both fireball and spec- 
tator, the separation between the two regions is clearer 
(smaller surface contacts) for lighter systems. 



normalized transverse momentum for each rapidity bin. 
It is evident that for the first two cases there is no scaling 
with respect to the system size, neither without nor with 
accounting for the system size via the (Ap /3 +A^ /3 ) term, 
while the in-plane average transverse momenta, propor- 
tional to vf, shows system size scaling (lower panel of 
Fig. 8) for the participant region. Somewhat expected, 
a deviation is present for the spectator part. We have 
observed a very similar feature for Z=3 particles, while 
for Z=l the scaling holds over all the forward rapidity 
domain. 
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FIG. 8. Integrated directed flow as a function of rapid- 
ity for Z=2 particles in the M4 centrality bin of collisions 
Au+Au, Xe+Csl and Ni+Ni at 250T MeV. Upper panel: vi 
values, middle panel: vi scaled by the term (A 1 ^ 3 + A^ 3 ), 
lower panel: scaled values, v" (see text). The lines are joining 
the symbols to guide the eye. 

Figure 8 shows the system dependence of the directed 
flow for Z=2 particles in the M4 centrality bin of colli- 
sions Au+Au, Xe+Csl and Ni+Ni at 250,4 MeV. Plotted 
are the integrated v\ values as a function of rapidity for 
three cases: i) as such (upper panel); ii) scaled with the 



term (Ap + A^ a ), which is proportional to the sum of 
radii of projectile and target (middle panel); iii) scaled as: 
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= vi(p[°^)/(Ap J + A^°), where {pt'') is the average 
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FIG. 9. Differential flow for three systems at 250,4 MeV, 
M4 centrality bin, for Z=2 particles in three windows of ra- 
pidity. The lines are polynomial fits to guide the eye. The 
arrows mark the values of the average p^for the correspond- 
ing system. 



The (A 1 / 3 



1 /3 

Arp ) scaling has been proposed by Lang 
et al. [34], who, within a BUU model, have found a linear 
dependence of the transverse pressure (leading to trans- 
verse momentum transfer) with the reaction time (pas- 
sage time). Westfall et al. [14] have related the A~ x l z de- 
pendence of Ebai to a competition between the attractive 
mean field (associated with the surface, so scaling with 
A 2 / 3 ) and the repulsive nucleon-nucleon interaction (scal- 
ing as A) . This competition between the two components 
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may be the origin of the quoted scaling of the transverse 
pressure with (Ap 3 + A}/ 3 ). Earlier studies devoted to 
the slope of (p x ) —y distributions (which translates into a 
flow angle) have experimentally confirmed such a scaling 
[18,22]. We note that, in an ideal hydrodynamics the flow 
angle is a pure geometric quantity and does not depend 
on the system size [56]. The system size dependences 
presented above may be an interesting effect of the nu- 
clear forces and/or a consequence of the non-equilibrium 
nature of the heavy-ion collisions. 

In Fig. 9 we present the differential flow for the three 
systems at 250A MeV, M4 centrality bin, for Z—2 parti- 
cles in three windows of rapidity (the three panels). The 
arrows mark the values of the average pj^for the corre- 
sponding centrality bin, according to the symbols. For 
this incident energy of 250,4 MeV the value of the projec- 
tile momentum in the cm. system is 342 MeV/c per nu- 
cleon. As seen already in Fig. 7, for all rapidity windows 
the (p^ ) depend on the system size (again breaking the 
scaling expected from hydrodynamics), suggesting an in- 
crease of the compression and expansion with the system 
size. 

Similar data as those presented in Fig. 8 and 9 are pre- 
sented in the Appendix for Z=l particles at 250A MeV 
and for Z=l and Z=2 particles at 400,4 MeV (Fig. 25 to 
Fig. 30). 

IV. MODEL COMPARISON 

The IQMD transport model [33,41] is widely used for 
interpreting the data in our energy domain [33,16,20,39]. 
We use two different parametrizations of the EoS, a 
hard EoS (compressibility K= 380 MeV) and a soft EoS 
(K= 200 MeV), both with MDI, labeled HM and SM, re- 
spectively and without MDI - H and S, respectively. We 
use the free nucleon-nucleon cross section, cr£^ ee for all 
cases, but for the energy of 90^4 MeV in addition we con- 
sider the case of a nn = 0.8c^ ee . The events produced 
by the model are filtered by the experimental filter and 
analyzed in a similar way as the experimental data. This 
comprises the same recipe for the centrality selection and 
the same way of reaction plane reconstruction and cor- 
rection. However, for the energy of 90A MeV, due to a 
weak flow signal (see below) we prefer to use the true 
reaction plane for the model calculations. In what con- 
cerns the reaction plane resolution at 400A MeV, it is in 
the model very similar compared to data. For instance, 
for M4 centrality bin, for IQMD SM the correction fac- 
tors are 1.03, 1.05 and 1.24 for Au, Xe and Ni systems, 
respectively. 



A. What to compare 

A known problem of the IQMD model (and QMD mod- 
els in general) is that of much lower yields of compos- 



ite fragments compared to data [57]. For instance, for 
Au+Au at incident energy of 400yl MeV, M4 central- 
ity, integrated Z=2 yields relative to Z=l are 1/4.7 for 
the experimental data while IQMD predicts 1/25 for HM 
and 1/15 for SM (these ratios do not depend on including 
MDI or not). 

As cluster formation and flow are intimately related, 
one cannot simply neglect this dramatic discrepancy. It 
is difficult to assess whether this strong dependence of 
fragment production on EoS parametrization is a gen- 
uine physical effect or is only particular to IQMD model. 
We note that most of the present models used in heavy- 
ion collisions at our energies involve a rather simple phase 
space coalescence mechanism to produce composite par- 
ticles [43,49] (IQMD uses the coalescence in coordinate 
space only [33]). Efforts to identify the fragments early 
in the collision might contribute to clarify this aspect 
[58] . Promising theoretical candidates to accomplish the 
task of realistic fragment formation could be AMD-type 
models [37]. 
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FIG. 10. Comparison of integrated v\ values as a function 
of rapidity for protons, neutrons, all particles weighted by 
mass (A) and all particles weighted by charge (Z) for IQMD 
HM events, for Au+Au at 400,4 MeV, M4 centrality. 

To partially overcome the problem of fragment produc- 
tion in the models, one can perform the comparison tak- 
ing into account all charged particles weighted by charge 
Z (so called proton-likes) [15,20]. However, this type of 
comparison could be biased, as the neutrons bound in 
the composite fragments may contribute differently for 
the calculations compared to data. We have investigated 
various possibilities by using IQMD events. The results 
are presented in Fig. 10, where we compare integrated v± 
values as a function of rapidity for protons, neutrons, all 
particles weighted by mass (A) and all particles weighted 
by charge (Z) for Au+Au at 400,4 MeV, M4 centrality 
bin. 
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FIG. 11. As Fig. 10, but for differential flow, in the rapidity 
window jp ' =0.7-0.9. The statistical errors are plotted for 
protons. 



In Fig. 11 we show the same comparison in case of dif- 
ferential flow. The neutrons exhibit the same flow as the 
protons, both for integrated and for differential values. 
As a consequence, within the model, in both cases the 
charge- weighted values are identical to the mass- weighted 
ones. However, as fragments heavier than Z=2 are ex- 
tremely few in the model, this result may be somewhat 
biased. In the following we are comparing data and 
model both for selected particle types and for proton- 
likes. 




FIG. 12. Integrated v\ values as a function of rapidity for 
Z=l and Z—2 particles, for the incident energy of 90^4 MeV. 
The data points (dots) are compared to IQMD calculations 
for two EoS (lines) parametrizations. The line labeled HM.8 
corresponds to HM case, using a nn — 0.8u^ ee . For the HM 
case the statistical errors of the model are plotted. 



B. Integrated values 

We start our comparison of the data with the IQMD 
model at the incident energy of 90A MeV, for the M4 
centrality bin. 

In Fig. 12 we show the rapidity dependence of the inte- 
grated v\ values for particles with Z=\ (upper panel) and 
Z=2 (lower panel). For both species the measured val- 
ues are compared to the IQMD calculations for HM and 
SM parametrizations. For the HM case an additional set 
of calculations has been performed using <j nn — 0.8er T {^ ee 
(labeled HM.8 in Fig. 12). For the HM case the statisti- 
cal errors are plotted. For the other cases the errors are 
comparable. For data the errors are in most cases smaller 
than the dimension of the points. The calculated values 
of the directed flow depend both on the parametrized 
EoS and, more pronounced, on a nn . This dependence 
is apparently of different magnitude for Z=l and Z=2 
particles. 

For the model calculations there is a coexistence of at- 
tractive (negative V\ values) and repulsive (positive vi) 
flow, manifested as a function of rapidity (we shall call 



this dual flow). This coexistence is different for the two 
particle species. We noticed that the above character- 
istics of the model calculations depend on centrality as 
well, both the magnitude of the dual flow and the parti- 
cle dependence being enhanced for more peripheral col- 
lisions. The model features are clearly not supported by 
the data, which show a monotonic repulsive flow over all 
the rapidity domain, both for Z=l and Z=2 particles, 
as seen in Fig. 12. For the experimental data, for the 
centrality M4 studied here, the reaction plane correction 
factor is 1.54. 

A two-component flow was observed earlier in QMD 
calculations of semi-peripheral Ca+Ca collisions at 
350 A MeV [39]. That study pointed out its high sen- 
sitivity to MDI. But, unless the discrepancy between the 
calculations and the measured data is resolved, any con- 
clusion on the sensitivity of the directed flow on the EoS, 
a nn or MDI is meaningless for energies around E^ai- It 
is not clear for the moment whether the particle depen- 
dence of the dual flow is not an artifact of the treatment 
of composite particles in the model. We note that mea- 
surements of Eb a i for different particle types [14,23] did 
not reveal, so far, any dual flow. Calculations with a 
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FIG. 13. Integrated V\ values as a function of rapidity, for 
Z=l particles, for three systems at 400A MeV, centrality bin 
M4. The data points (dots) are compared to IQMD calcula- 
tions (lines). 

BUU model [44] found a dual flow only in a (p t ) differ- 
ential way, but otherwise monotonic behavior of (px) — V 
distributions. Recent experimental investigations of flow 
in light systems at Eb a i pointed out interesting aspects 
of flow of light isotopes and heavy fragments, but again 
the balance energy was found not to depend on particle 
type [59]. 

In Fig. 13 and Fig. 14 we show the comparison of data 
and model calculations for Z=l and Z=2 particles, re- 
spectively. The three studied systems are considered, for 
the centrality bin M4. The integrated v\ values show 
sensitivity to both EoS and to MDI. As expected, the 
MDI influence the flow essentially in the vicinity of the 
projectile spectator (j/°) > 0.8). This effect is more pro- 
nounced the lighter the system. All these sensitivities 
are enhanced for Z=2 particles. For Au+Au and Xe+Csl 
systems, the SM parametrization is reproducing the data 
very well, for both Z=\ and Z—2 particles. This may be 
the result of a similar balance of thermal and collective 
contributions in the model compared to the data. In fact, 
the phase space populations of Z—\ and Z—2 particles 
are similar for model and data. 
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FIG. 14. As Fig. 13, but for Z=2 particles. For the HM 
case the statistical errors of the model are plotted. 

In case of Ni+Ni system the dependence on EoS is al- 
ready negligible, but obviously the model underestimates 
the flow. One parameter of the model, the Gaussian 
width L, which is the phase space extension of the wave 
packet of the particle (and acting as an effective interac- 
tion range) has been found to influence the directed flow 
considerably [41]. A decrease of L for lighter systems has 
been advocated with the argument of maximum stability 
of nucleonic density profiles [41]. As no clear prescription 
exists for handling the value of L, we prefer to use a con- 
stant value of L=8.66 fm 2 throughout the present work. 
A smaller L would lead to an increase of the directed flow 
[41] and may cure the discrepancy that we observe for the 
Ni+Ni system, but will affect unfavorably the compari- 
son in case of Xe+Csl system. These effects may reflect 
the importance played (via the interaction range) by the 
surface. A complete understanding of this aspect is a 
necessary step towards establishing the bulk properties 
of the nuclear matter created in heavy-ion collisions. 

We note that comparisons of integrated directed flow 
for Au+Au system using QMD-type models favored 
mostly a soft EoS [15,54] but a hard EoS was also found 
to explain another set of experimental data [16]. 

In Fig. 15 we present the comparison of the measured 
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FIG. 15. Integrated V\ values as a function of rapidity for 
all particles weighted by Z, for the centrality bins M4 (up- 
per panel) and M3 (lower panel), for the incident energy of 
400,4 MeV. The data points (dots) are compared to IQMD 
calculations (lines). 



integrated v% values to IQMD calculations for Au+Au at 
the incident energy of 400,4 MeV, taking into account all 
charged particles weighted by charge Z. The centrality 
bins M4 and M3 are studied. In this case the sensitivity 
to EoS is reduced, as a consequence of a balance between 
magnitude of flow and yield of composite particles in the 
model: hard EoS produces more flow, but less particles 
with Z >1, while for soft EoS it is opposite. This be- 
havior strongly underlines once more the necessity that 
theoretical models appropriately describe the yields of 
composite particles. The conclusion on EoS is this time 
less evident, but the parametrizations without MDI are 
ruled out once again, on the basis of their departure from 
the data in the region of spectator rapidity. As expected, 
this effect is more pronounced for the more peripheral 
centrality bin M3. 

Despite the good agreement seen at the beam energy 
of 400,4 MeV, we found that in the IQMD model the 
decrease of flow towards lower incident energies is much 
faster than for data, leading to larger theoretical Eb a i 
compared to data (and to the behavior seen in Fig. 12). 
This may be a result of deficiencies in incorporating MDI 
and in the treatment of fragment production. The Pauli 
blocking may play a role too. In addition, it has been 



pointed out that the shape of the flow excitation func- 
tion is drastically influenced by the method of imposing 
constraints on the Fermi momenta [41]. The features 
of the model calculations presented above for 400,4 MeV 
show the danger of deriving EoS-related conclusions from 
rapidity- integrated flow values (like p* r ) unless a de- 
tailed description of the data is first achieved in a dif- 
ferential way. As realized early on [27], a soft EoS with 
MDI is producing similar magnitude in p^ r as a hard 
EoS without MDI. 



C. Differential flow 

We restrict our model comparison of the differential 
flow to the incident energy of 400,4 MeV and M4 cen- 
trality bin. Data for all three systems investigated so far 
are compared to the model calculations. 
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p/0) 

FIG. 16. Differential flow for particles with Z=l (upper 
row) and Z—2 (lower row) for Au+Au collisions at incident 
energy of 400 A MeV, M4 centrality, for two windows in ra- 
pidity (columns). Experimental data are represented by dots 
and the model calculations are the lines. For the HM case the 
statistical errors of the model are plotted. 

In Fig. 16 the measured differential directed flow for 
Au+Au collisions at incident energy of 400,4 MeV, M4 
centrality, is compared to the IQMD results for all the 
four parametrizations used above. Particles with Z—l 
(upper row) and Z=2 (lower row) for two windows in 
rapidity are used for the comparison. For both particle 
species there is a clear sensitivity of DDF on the EoS. As 
for the case of the integral flow, the SM parametrization 
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reproduces the experimental data quite well. Apparently 
the model calculations deviate from the data at high p t 
in case of Z=l particles, while the corresponding Z=2 
particles are well explained. This deviation is more pro- 
nounced for larger rapidities. We have found earlier [60] 
that a BUU model does not explain the DDF of protons 
in the spectator region at higher energies. The shape 
of the Z=l DDF distributions are in case of the IQMD 
model strikingly similar to the ones of Z=2, while for the 
data there are subtle differences between the two parti- 
cle species (at this energy of 400^4 MeV as well as down 
to 90 A MeV [48]). The model features may be the re- 
sult that the nucleons (dominating the Z—l sample) in 
the models are all "primordial" , which does not account 
for the sequential decays of heavier fragments. The dy- 
namics of the expansion and fragment formation may be 
responsible for the differences, too. 

TABLE III. Average normalized transverse momentum 
(Pt ) for particles with Z=l and Z=l in Au+Au collisions 
at 400A MeV, M4 centrality bin. Data and model values are 
compared for two rapidity windows. For data, the number in 
parenthesis represents the error on the last digit. 
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FIG. 17. Comparison of data and model differential flow 
for Z—l particles, for the incident energy of 400A MeV, M4 
centrality bin, for the systems Xe+Csl (upper row) and Ni+Ni 
(lower row) for two windows in rapidity (columns). 



In Table III we compare the experimental values of 
the average normalized transverse momentum with the 
values from IQMD, for HM and SM cases. Particles with 
Z=\ and Z—2 for the two windows in rapidity studied in 
Fig. 16 are compared. The data values have a systematic 
error represented by the number in parenthesis as the 
error on the last digit. The model reproduces reasonably 
well the average transverse momenta for Z=l particles, 
while it underestimates them for Z—2, for both windows 
of rapidity. 

We mention that, recently, our experimental differ- 
ential directed flow in Au+Au [48] was nicely repro- 
duced by a BUU model wich includes an improved Dirac- 
Brueckner formalism [49]. In this case, for the densities 
expected at 400 A MeV, the EoS is soft, which is in agree- 
ment to our results. 

In Fig. 17 we show the measured DDF for Xe+Csl and 
Ni+Ni systems at 400A MeV, M4 centrality, in compar- 
ison to IQMD results, for Z=l particles in two windows 
of rapidity. In case of Xe+Csl system the model calcula- 
tions are at the same level of agreement with data as in 
case of Au+Au: the SM parametrization is reproducing 
the data, with clear deviations at high momenta. For the 
Ni+Ni case even the HM parametrization underpredicts 
the measured data. Most notably, as obvious particu- 
larly for Ni+Ni, the MDI have effects predominantly at 
low pt, contrary to earlier BUU predictions (performed 
for the asymmetric system Ar+Pb) [36]. 



In Fig. 18 we show the model comparison of the differ- 
ential flow for particles weighted by Z for Au+Au colli- 
sions at incident energy of 400A MeV, M4 centrality, for 
two windows in rapidity. As in case of integrated values, 
as a result of different relative contribution of particles 
heavier than Z=l, the sensitivity to EoS is reduced for 
this type of comparison. 

V. SUMMARY AND CONCLUSIONS 

We have presented experimental results on directed 
flow in Au+Au, Xe+Csl and Ni+Ni collisions at incident 
energies from 90 to 400A MeV. General features of the 
directed flow have been investigated using experimental 
data, particularly the centrality and the system depen- 
dence. We have studied the rapidity dependence of the 
first Fourier coefficient, V\, integrating over all transverse 
momentum range. A special emphasis has been put on 
the differential directed flow, namely the pt dependence 
of v\ . While for integrated values we presented a new way 
of looking at old (and generally known) dependences, the 
DDF results are reported for the first time for our energy 
domain, both for the centrality and for the system size 
dependence. We have devoted special care to the cor- 
rections of the experimental data. The influence of the 
finite granularity of the detector has been studied and 
corrected for. The high accuracy of the final results is 
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FIG. 18. Differential flow for particles weighted by Z for 
Au+Au collisions at incident energy of 400 A MeV, M4 cen- 
trality, for two windows in rapidity. Experimental data are 
represented by dots. The model calculations are the lines. 



tion of system size, and iii) fragment production. As a 
consequence, none of the IQMD parametrizations stud- 
ied here is able to consistently explain the whole set of 
experimental data. 

The importance of spectators acting as clocks for the 
expansion is one particular argument to study collective 
flow in semi-central collisions at energies from a few hun- 
dred MeV to a few GeV per nucleon [4] . We have demon- 
strated that high precision experimental data allows us to 
study the many facets of the heavy-ion collisions. Other 
observables, like 1)2, should receive a comparable (and 
simultaneous) attention too. Whether the nuclear equa- 
tion of state can be extracted from such studies depends 
ultimately on the ability of any type of microscopic trans- 
port model to reproduce the measured features. 
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based as well on a good reaction plane resolution achieved 
with the full coverage of the FOPI detector. 

We have compared the experimental data with IQMD 
transport model calculations, for both integral and differ- 
ential v\ values. This comparison, performed for all the 
three studied systems, shows a clear sensitivity of the 
directed flow on the EoS parametrization in the model, 
especially in case of particle-selected comparison. In this 
case, for both integrated and differential directed flow 
at the incident energy of 400^4 MeV, we conclude that 
a soft EoS with MDI is the only parametrization in the 
model that reproduces the data for Au and Xe systems. 
A clear discrepancy is seen for Ni system, which needs 
to be addressed separately. It may reflect the increasing 
importance played by the nuclear surface for lighter sys- 
tems. We consider our present results as a case study on 
the sensitivities in determination of EoS and MDI from 
directed flow comparisons. We emphasized the necessity 
of the present kind of differential comparison prior to 
more global quantities. We have shown that the combi- 
nation of rapidity and transverse momentum analysis of 
(differential) directed flow can impose constraints on the 
model. We also pointed out some difficulties of the model 
to reproduce the measured data concerning: i) flow at low 
energy (we considered here 90 A MeV) , ii) flow as a func- 



13 



APPENDIX 



In the following we present additional results that com- 
plement those included in the body of the paper. As the 
trends are similar to those already discussed in section 
III, we present here only figures. 




FIG. 20. Centrality dependence of the integrated directed 
flow as a function of rapidity for Au+Au at 250A MeV for 
Z=2 particles. 
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FIG. 21. Centrality dependence of the integrated directed 
flow as a function of rapidity for Au+Au at 400^4 MeV for 
Z=\ particles. 




p/o) 

FIG. 22. Differential flow for three centrality bins, in three 
rapidity windows, for Z=l particles for collisions Au+Au at 
25CL4 MeV. The lines are polynomial fits to guide the eye. 
The arrows mark the values of the average p\ for the corre- 
sponding centrality bin. 
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FIG. 23. Differential flow for three centrality bins, in three 
rapidity windows, for Z—2 particles for collisions Au+Au at 
25(L4 MeV. 



FIG. 24. Differential flow for three centrality bins, in three 
rapidity windows, for Z=l particles for collisions Au+Au at 
40CL4 MeV. 
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FIG. 25. Integrated directed flow as a function of rapid- 
ity for Z=\ particles in the M4 centrality bin of collisions 
Au+Au, Xe+Csl and Ni+Ni at 250,4 MeV. Upper panel: vi 
values, middle panel: vi scaled by the term (A^ 3 + Aj! 3 ), 
lower panel: scaled values, vf = vi (p[°^ ) /(Ap 3 + A 1 ^ 3 ). 




FIG. 26. Integrated directed flow as a function of rapid- 
ity for Z=\ particles in the M4 centrality bin of collisions 
Au+Au, Xc+Csl and Ni+Ni at 40(L4 MeV. 
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FIG. 27. Integrated directed flow as a function of rapid- 
ity for Z=2 particles in the M4 centrality bin of collisions 
Au+Au, Xe+Csl and Ni+Ni at 40(L4 MeV. 



FIG. 28. Differential flow for three systems at 250,4 MeV, 
M4 centrality bin, for Z=l particles in three windows of ra- 
pidity. The lines are polynomial fits to guide the eye. The 
arrows mark the values of the average p)°'ior the correspond- 
ing system. 
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FIG. 29. Differential flow for three systems at 400,4 MeV, 
M4 centrality bin, for Z—l particles in three windows of ra- 
pidity. 
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FIG. 30. Differential flow for three systems at 400,4 MeV, 
M4 centrality bin, for Z=2 particles in three windows of ra- 
pidity. 
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